Bose-Einstein condensate in gases of rare atomic species 
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We report on the successful extension of production of Bose-Einstein Condensate (BEC) to rare 
species. Despite its low natural abundance of 0.13%, ^^*Yb is directly evaporatively cooled down 
to BEC. Our successful demonstration encourages attempts to obtain quantum gases of radioactive 
atoms, which extends the possibility of quantum many-body physics and precision measurement. 
Moreover, a stable binary mixture of ^''^Yb BEC and ^^*Yb BEC is successfully formed. 
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Successful realization of Bose-Einstein condensates 
(BEC) and Fermi degenerate gases in atomic vapors 
has been providing intriguing possibilities in many re- 
search fields such as the study of quantum many-body 
system, precision measurement, and quantum compu- 
tation(see Ref.lJ for review). So far atomic species 
with relatively high natural abundance have been Bose- 
Einstein condensed, such as in alkali-atoms®^'®^Rb[2, yjj 
39^4iK[4, 5], "Li 6], 23Nala] and ^^^CstS], and two-electron 
atoms such as ^Hei, [if, ^'^'^'™'™Yh^MM: *°Ca[l|, 
and 84^86,88gj.[J^|g^^ and also ^H^ and 52Cr[23l. Ex- 
tending atomic species beyond these relatively abundant 
atomic species is an important step for a future investi- 
gation. First, extension to rare species benefits a richer 
variety of quantum degenerate gases and the mixed gases 
in the study of quantum many-body system. In ad- 
dition, radioactive atoms such as francium[2l|, Um ^^.d 
radium |23| are good candidates for tests of the fundamen- 
tal symmetries by searching for the electric dipole mo- 
ment and precision measurement of atomic parity non- 
conservation. In such application, a Mott-insulator state 
in an optical lattice, realized by making use of quantum 
gases, offers important advantage of high atomic density 
with avoided collision by localizing each atom to each 
lattice site. 

Obviously, the limited abundance of atoms is a great 
obstacle to reach quantum degeneracy [24J. However, in 
principle, there is no fundamental difficulty to achieve 
quantum degeneracy of rare species by evaporative cool- 
ing if the ratio of elastic to inelastic scattering cross sec- 
tion is high. In this Rapid Communication, we report 
on the successful production of BEC of rare species of 
^^*Yb, which has the lowest natural abundance of only 
0.13% among the seven stable isotopes of Yb without the 
use of enriched source. The abundance of 0.13% is also 
the lowest among the already produced BEC 2J| . The 
atom clouds can be directly evaporatively cooled down 
to BEC owing to a favorable s-wave scattering length of 
^^*Yb and a high vacuum of our system. The lowest 
natural abundant species so far was ^"^Sr, which has a 
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natural abundance of 0.56%. The difference from ^''Sr 
BEC is the small initially-captured atom number. While 
10^ atoms were magnetically trapped in ^^Sr, only 2 x 10^ 
atoms are captured in a magneto-optical trap (MOT) for 
^^^Yb. This value is, for example, not too far from the 
atom number of radioactive francium MOT[25|. There- 
fore, our successful demonstration encourages attempts 
to obtain quantum gases of radioactive atoms [26j. A 
straightforward implication of our result is the promis- 
ing possibility of realizing quantum degenerate gases of 
radioactive Yb isotopes for precision measurement '27l|. 

In addition, we report on the attainment of dual- 
species condensates of ^^^Yb and ^^^Yb. This is the first 
stable binary condensate for a two-electron system. Since 
two species BEC has been reported only for ^^Rb and 
s^Rb[28l, "^K and 8^Rb[29i, and quite recently «^Rb and 
^■^■^Cs mixtures [30], the demonstrated dual condensate of 
^^*Yb and ^''^Yb in fact extends the variety of the quan- 
tum degenerate mixture, where quantum turbulences |31j 
and topological excitations [33 can be studied. 

In order to produce ^^®Yb BEC, we have revised our 
experimental system from the previously reported ones 
[Il| - Kl3l . l33j . The major improvement is the high vac- 
uum of l~3xl0~^^ Torr compatible with the highly 
coUimated atomic beam. The experimental procedure 
to produce ^^^Yb BEC follows our previous scheme (lll- 
\13i . |33| . |34| . First, an ^^®Yb atomic beam from an oven 
containing Yb metal with natural abundance is deceler- 
ated by a Zeeman slower operating on the ^Sq -h- ^Pi 
transition with a laser power of 250 mW. The deceler- 
ated atoms are collected by a MOT operating on the 
narrow-line inter-combination transition of ^Sq O '^Fi. 
Our high vacuum system enables us to typically collect 
2x 10^ 168 YI3 atoms at a temperature of 20 /iK during the 
100 s loading. This should be compared with the case of 
the most abundant isotope of ^^'*Yb where it takes less 
than 0.5 s to collect 2 x 10® atoms. The loading time of 
about 100 s would be limited by a very little branching 
from the "^Pi state to the '^Pq state throu gh t he magnetic 
dipole transition at a rate of '^ 1/161 s~"'^[33 and also by 
background gas collisions. 

Then we transfer the atoms from the MOT to a Far- 
Off Resonant Trap (FORT) operating at 532 nm to per- 



form efRcient evaporative cooling. Efficient transfer is 
performed by compressing the magnetic field gradient of 
the MOT and subsequently reducing the MOT intensity 
and the detuning of MOT laser beams to cool the atoms. 
Typically 1x10^ ^^*Yb atoms are successfully trapped 
into the FORT. 

The evaporative cooling is performed in a crossed 
FORT, where horizontally and vertically propagating 
FORT beams intersect at their foci. The horizontally 
propagating FORT beam is elliptically shaped with a 
beam waist of 30 fira along the horizontal direction and 
14 /^m along the vertical direction. The initial trap depth 
is approximately 750 /xK at the horizontal FORT laser 
power of 14 W. The first stage of evaporative cooling is 
performed by lowing the laser intensity of the horizontal 
FORT beam, during which the atoms are transfered to 
the crossed FORT region. The shape of the vertically 
propagating FORT beam is also elliptical with a beam 
waist of 127 /im along the axial direction of the horizon- 
tal FORT beam and 40 /im along the radial. This results 
in a good overlap of the two FORT beams with an in- 
creased trap volume. The potential depth by the vertical 
FORT beam is 40 ^K. At 40 fiK potential depth of the 
crossed FORT, 5 x 10^ atoms at 6 /iK are trapped in the 
crossed region. 

The atoms in the crossed FORT region are further 
cooled down toward the quantum degenerate regime 
at the final stage of the evaporative cooling. The s- 
wave scattering length of ^^^Yb is determined to be 
0-168 = 13.33(18) nm from two-photon photoassociation 
spectroscopy [36|. The evaporative cooling should work 
well at this scattering length. To suppress unfavorable 
inelastic three-body loss, we not only reduce the hori- 
zontal FORT beam intensity, but also the vertical FORT 
beam intensity during the evaporation to reduce the cen- 
tral peak density. The peak density above the BEC phase 
transition is typically 2^4x10^'^ cm~^. 



The phase transition from a thermal gas to a BEC 
is observed as the appearance of bimodal distribution 
after Time-of-Flight (TOF). Figure 1 shows absorption 
image taken after TOF of 14 ms for different evapo- 
ration stage near BEC phase transition. At a higher 
temperature, the density distribution shows an isotropic 
Gaussian distribution, manifesting that the velocity of 
the gas has a Maxwell-Boltzmann distribution. How- 
ever, after further evaporative cooling, we observe the 
deviation from the Gaussian distribution and an addi- 
tional low velocity component, representing a BEC com- 
ponent. At this stage, the peak density of the thermal 
gas is 3.5 X 10^^ cm~'^. The trap frequencies are 2 tt 
X (240, 31 165) Hz. The transition temperature is eval- 
uated to be Tc = 135 nK. This value is slightly lower 
than the transition temperature for non-interacting BEC, 
Tc = Uuj/kB X (iV/C(3))i/3 ^ 150 nK, where w is the ge- 
ometric mean of trap frequency, C is the Riemann zeta 
function, and N is the atom number. This slight de- 
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FIG. 1: (Color online) Phase transition to BEC during the 
evaporative cooling, (a) Absorption images at different evap- 
oration stage are shown with TOF of 14 ms. Each image is 
290 lira per side, (b) The corresponding integrated column 
densities along the horizontal direction are shown with laser 
powers of horizontal FORT beam. The BEC transition is ob- 
served as an appearance of bimodal distribution (middle). At 
the end of the evaporation, a almost pure condensate of 1.0 x 
10^ is produced (top). The absorption images are taken from 
an angle of 45 degree against the horizontal FORT axis in the 
horizontal plane. 



viation is expected from the interaction and finite size 
effects [131. 

After further evaporation by lowing the trap depth, the 
condensate fraction increases. After evaporative cooling 
for 9.3 s, we finally attain an almost pure condensate 
up to 1.0 X 10^ atoms. The trap frequencies are 2 tt 
X (154, 25, 141) Hz. The peak density is 1.0 xlO^^ 
cm^'^. The validity of Thomas- Fermi approximation is 
confirmed by the relation Naies/cho — 1-4 x 10^ ^ 1, 
where Oho = yhJijmJj) is a harmonic oscillation length 
and TO is the atomic mass. Figure 2 shows radii of the 
expanding BEC during the TOF. Anisotropic expansion 
of condensate can be clearly seen. The solid lines in Fig.l 



are the fits using the Castin-Dum model with the chemi- 
cal potential as the only free parameter. The correspond- 
ing chemical potential is 64 nK. It is noted that, with this 
achievement, all the isotopes of Yb are now cooled down 
to quantum degenerate regime except for the bosonic iso- 
tope of ^^^Yb, where a large number of condensate cannot 
grow due to its large negative scattering length [36]. 

At the final trap depth, a lifetime of BEC is evalu- 
ated to be about 3 s, which is mainly attributed to the 
three-body loss of atoms. From this, we can determine 
the upper limit of three-body loss rate of ^^^Yb. We fit 
the remaining atom number as a function of holding time 
in the trap by a rate equation assuming only for three- 
body loss. The maximum three-body coefficient K^ is 
evaluated to be 8.6(1.5) x 10~^* cm^/s. This is more 
than one order of rnagnitude larger than the value re- 
ported for ^^"'Ybjllllljl. However it is still smaller than 
the expected maximal rate of Ks^max — 210haiGs^/m ~ 

2.3X 10"27(,j^6/g^^J^ 
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FIG. 2: (Color online) Expansion of condensate after release 
from the trap. The radii of the condensate show asymmetry 
between the horizontal (circle) and vertical (square) direction. 
The solid lines are the fits using Castin-Dum model. The error 
bars are standard error. 



factor of 230, we load ^^^Yb for 100 s and then load ^''^Yb 
for 0.5 s by subsequently turning on the Zeeman slower 
beam. Since the interspecies scattering length between 
the two species is as small as 0.13(18)nm, both species are 
almost independently evaporatively cooled to quantum 
degenerate regime. Figure 3 shows the absorption images 
of the dual condensate. The number of atoms in each con- 
densate is 0.9x10''. The advantage of working with Yb 
isotope mixture is the fact that trap parameters such as 
trap depth, trap frequency, and gravitational sag, are al- 
most the same for the different isotopes. This enables us 
to focus on the role of interaction. Such feature is already 
seen in the behavior of the expanding condensate in Fig. 
3, where ^^^Yb is expanding much faster than ^^^Yb. As 
the interspecies interaction is small, the difference in the 
size can be explained by the mean field energy of each 
condensate, which is converted into kinetic energy dur- 
ing the expansion. The expansion energy is conserved to 
2fi/7, which is proportional to (NaM)^^^ in the Thomas- 
Fermi approximation, where N is the number of conden- 
sate, ^ is the chemical potential and om is the scattering 
length with mass number M. The expansion velocity 
normalized by N^^^ of ^^®Yb condensate is evaluated to 
be 1.20(4) times faster than ^^^Yb from the TOF image 
of each condensate. This is in good agreement with the 
ratio of scattering length (ai68/ai74)^^^=l-19. The ratio 
of each number of atoms can be easily tuned by chang- 
ing the loading time to the MOT. We can so far prepare 
an almost pure Bose-Bose mixture with ^^^Yb up to 7 x 
10^. Although the s-wave interspecies scattering length 
is small, optical Feshbach resonance will enable tuning of 
interspecies interaction, which is demonstrated for BEC 
of i^^'YbfH and thermal gases of "^yb and i^^Ybfill. 
Study of quantum quench dynamics and non-equilibrium 
dynamics [431 in binary condensate system will also be ac- 
cessible. 



Our successful demonstration of BEC of rare species 
extends the possibility to create various quantum degen- 
erate mixtures. This is especially important for Yb quan- 
tum isotope mixture, since a stable BEC mixture has 
not been yet reported. Previously, we reported on the 
production of Bose-Bose mixture of ^^'^Yb-^^^Yb below 
Tc[1j1. However, due to the negative scattering length of 
^^^Yb, a condensate above the critical atom number (typ- 
ically a few hundred) gets unstable. Among the possible 
quantum degenerate Bose-Bose mixtures of Yb isotope, 
i68Yb-i74Yb and i^Syb-i^Oyb ^re the only candidate to 
realize a stable BEC mixture. Here, we report on the 
production of stable ^®^Yb-^^^Yb Bose-Bose mixture. 

The procedure to produce this mixture is basically the 
same as that describe above. To collect two isotopes in 
the MOT, we apply two-color MOT beam. Due to the 
large difference of the atomic beam fiux, originates from 
the large difference in the natural abundance of about a 
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FIG. 3: (Color online) A stable Bose-Bose mixture of ^'^^Yb- 
^■^''Yb. Absorption images of (a) almost pure "*Yb BEC of 
9x10^ and (b) almost pure "*Yb BEC of 9x10^. Both images 
are taken after TOF time of 18 ms in a separate experimental 
run. The larger cloud size of ^^*Yb compared with ^^^Yb 
manifests larger mean field energy of the condensate of ^^* Yb 
compared with ^^^Yb. 



In conclusion, we have produced a BEC of rare species 
of ^^*Yb containing 1.0 x 10'* atoms. With this achieve- 
ment all the possible isotopes of Yb have been cooled to 
quantum degeneracy. Despite its low natural abundance 
of 0.13%, ^^^Yb is efficiently collected and cooled down to 
BEC. Our demonstration encourages researchers that a 
radioactive isotope is a realistic candidate of BEC. More- 
over, is^Yb BEC is produced together with ^'^''Yb BEC 
as a stable binary BEC. This mixture, which has very 
small interspecies scattering length, will be a prime can- 



didate for studying a binary BEC, including the possible 
use of optical Feshbach resonances. 
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